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1 2D and 3D Heat Equation 

Ref: Myint-U & Debnath §2.3 - 2.5 
[Nov 2, 2006] 

Consider an arbitrary 3D subregion V of M 3 (V C IR 3 ), with temperature u (x, t) 
defined at all points x = (x, y, z) G V . We generalize the ideas of 1-D heat flux to 
find an equation governing u. The heat energy in the subregion is defined as 



heat energy = J J cpudV 

Recall that conservation of energy implies 

rate of change heat energy into V from heat energy generated 
of heat energy boundaries per unit time in solid per unit time 

We desire the heat flux through the boundary S of the subregion V, which is 
the normal component of the heat flux vector 0, • n, where n is the outward unit 
normal at the boundary S. Hats on vectors denote a unit vector, |n| = 1 (length 1). 
If the heat flux vector is directed inward, then <fi ■ n < and the outward flow of 
heat is negative. To compute the total heat energy flowing across the boundaries, 
we sum • n over the entire closed surface S, denoted by a double integral J J s dS. 
Therefore, the conservation of energy principle becomes 

iIII v ^ dV = -II/ &dS+ III v QdV (1) 
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1.1 Divergence Theorem (a.k.a. Gauss's Theorem) 

For any volume V with closed smooth surface S, 




V- AdV = / / A-hdS 

v J Js 



where A is any function that is smooth (i.e. continuously differentiable) for x e V. 
Note that 

/ d d d \ „ d „ d „ d 
ydx'dy'dz) x dx y dy z dz 

where e x , e y , e z are the unit vectors in the x, y, z directions, respectively. The 
divergence of a vector valued function F = (F x , F y , F z ) is 

dF x , dF y , dF z 



V-F = — - + — A + 



dx dy dz 

The Laplacian of a scalar function F is 

9 2 F 9 2 F 9 2 F 
V 2 F = 1 1 . 

dx 2 dy 2 dz 2 

Applying the Divergence theorem to (1) gives 

d 
dt 




cpudV = - V-(j)dV+ / / / QdV 

\ J J Jv J J Jv 



Since V is independent of time, the integrals can be combined as 

du \ 
cp— +V-0-QJ dV = 




v 



Since V is an arbitrary subregion of M 3 and the integrand is assumed continuous, the 
integrand must be everywhere zero, 

du 

cp— + V • - Q = (2) 

1.2 Fourier's Law of Heat Conduction 

The 3D generalization of Fourier's Law of Heat Conduction is 

= -K Vu (3) 
where K is called the thermal diffusivity. Substituting (3) into (2) gives 



where k = K / (cp). This is the 3D Heat Equation. Normalizing as for the ID case, 

X K 
X=y, t=-t, 

Eq. (4) becomes (dropping tildes) the non-dimensional Heat Equation, 
where q = l 2 Q/ (ncp) = l 2 Q/K . 



2 2D and 3D Wave equation 

The ID wave equation can be generalized to a 2D or 3D wave equation, in scaled 
coordinates, 

u tt = V 2 m (6) 

This models vibrations on a 2D membrane, reflection and refraction of electromagnetic 
(light) and acoustic (sound) waves in air, fluid, or other medium. 

3 Separation of variables in 2D and 3D 

Ref: Guenther & Lee §10.2, Myint-U & Debnath §4.10, 4.11 

We consider simple subregions DCK 3 . We assume the boundary conditions are 
zero, u = on dD, where dD denotes the closed surface of D (assumed smooth). 
The 3D Heat Problem is 

u t = V 2 m, x G D, t>0, 
«(x,t) = 0, xedD, (7) 
«(x,0) = /(x), xGD. 

The 3D wave problem is 

u tt = V 2 m, x e D, t > 0, 

u(x,i) = 0, xe&D, (8) 

u(x,0) = /(x), xeA 

u t (x, 0) = g (x) , x£D. 



We separate variables as 



u (x, i) = X (x) T (t) (9) 
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The 3D Heat Equation implies 

T' W^X 

— = ~y~ = -A = const (10) 

where A = const since the l.h.s. depends solely on t and the middle X"/X depends 
solely on x. The 3D wave equation becomes 

T" V 2 X 

— = ~y~ = -A = const (11) 

On the boundaries, 

X (x) = 0, x G 3D 
The Sturm-Liouville Problem for X (x) is 

V 2 X + AX = 0, xefl (12) 
X(x) = 0, x G 3D 

4 Solution for T (t) 

Suppose that the Sturm-Liouville problem (12) has eigen-solution X n (x) and eigen- 
value A n , where X n (x) is non-trivial. Then for the 3D Heat Problem, the problem 
for T(t) is, from (10), 

V 

Y = ~ X (13) 

with solution 

T n (t) = c n c- x ^ (14) 
and the corresponding solution to the PDE and BCs is 

u n (x, t) = X n (x) T n (t) = X n (x) c n e- Xnt . 

For the 3D Wave Problem, the problem for T (t) is, from (11), 

rpll 

Y = ~ x ( 15 ) 

with solution 

T n (t) = a n cos {\/~\nt ) + p n sin ^A/A^t j (16) 
and the corresponding normal mode is u n (x, t) = X n (x) T n (t). 
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5 Uniqueness of the 3D Heat Problem 



Ref: Guenther & Lee §10.3 

We now prove that the solution of the 3D Heat Problem 

u t = V 2 w, x G D 
u(x,t) = 0, xG&D 
«(x,0) = /(x), X6D 

is unique. Let u±, 112 be two solutions. Define v — u\ — u-i- Then v satisfies 

v t = V 2 v, x G D 
v (x, t) = 0, x G dD 
u(x,0) = 0, xeD 

Let 

V{t)= J J J v 2 dV > 

V (t) > since the integrand f 2 (x, t) > for all (x, t). Differentiating in time gives 

dV 




dt (t) = I I I 2vv t dV 



D 



Substituting for v t from the PDE yields 

dV 
~dt 

By result (26) derived below, 
dV 



2vV 2 vdV 




(t) = 2 I I vVv-MS-2 I I I \Vv\ 2 dV (17) 



<dD J J JD 

But on dD, v — 0, so that the first integral on the r.h.s. vanishes. Thus 

dV 
~dt 

Also, at t — 0, 




(/) = -2 / / / \Vv\ 2 dV < (18) 



v (o) = y y y w 2 ( X , 0)^=0 



Thus V (0) = 0, V (t) > and dV/dt < 0, i.e. V (t) is a non-negative, non-increasing 
function that starts at zero. Thus V (t) must be zero for all time t, so that v (x, t) 
must be identically zero throughout the volume D for all time, implying the two 
solutions are the same, u\ — U2- Thus the solution to the 3D heat problem is unique. 

For insulated BCs, Vi> = on dD, and hence vVv ■ n = on dD. Thus we can 
still derive Eq. (18) from (17), and the uniqueness proof still holds. Thus the 3D 
Heat Problem with Type II homogeneous BCs also has a unique solution. 
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6 Sturm-Liouville problem 



Ref: Guenther & Lee §10.2, Myint-U & Debnath §7.1 - 7.3 

Both the 3D Heat Equation and the 3D Wave Equation lead to the Sturm-Liouville 
problem 

V 2 X + AX = 0, xg£>, (19) 
X(x) = 0, x G 3D. 

6.1 Green's Formula and the Solvability Condition 

Ref: Guenther & Lee §8.3, Myint-U k Debnath §10.10 (Exercise 1) 

For the Type I BCs assumed here (u (x, t) = 0, for x G dD), we now show that 
all eigenvalues are positive. To do so, we need a result that combines some vector 
calculus with the Divergence Theorem. From vector calculus, for any scalar function 
G and vector valued function F, 

V • (GF) = GV • F + F • VG (20) 
Using the divergence theorem, 

J J (GF)-hdS = J J J V-{GF)dV (21) 
Substituting (20) into (21) gives 

J J {GF)-MS = J J J GV -FdV + J J J F-VGdV (22) 
Result (22) applied to G — Vi and F = Vi> 2 gives 

j J (viVv 2 ) ■ MS = J J J ( Vl V 2 v 2 + Vv 2 ■ Vui) dV (23) 
Result (22) applied to G — v 2 and F = Vi>i gives 

J J {v 2 Vv 1 )-hdS = J J J (v 2 V 2 V! + Vui • V^ 2 ) dV (24) 

Subtracting (23) and (24) gives Green's Formula (also known as Green's second iden- 
tity): 

J J {v l Vv 2 -v 2 Vv 1 )-MS = J J J ( Vl V 2 v 2 - v 2 VV) dV (25) 

This is also known as a Solvability Condition, since the values of V\ and v 2 on the 
boundary of D must be consistent with the values of V\ and v 2 on the interior of 
D. Result (25) holds for any smooth function v defined on a volume V with closed 
smooth surface S. 
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6.2 Positive, real eigenvalues (for Type I BCs) 

Ref: Myint-U & Debnath §7.2 

Choosing G = v and F = Vv, for some function v, in (22) gives 





vVv-MS =111 (vV 2 v + Vv-Vv)dV 

(vV 2 v + \Vv\ 2 ) dV. (26) 




v 



Result (26) holds for any smooth function v defined on a volume V with closed smooth 
surface S. 

We now apply result (26) to a solution X (x) of the Sturm-Liouville problem (19). 
Letting v = X (x), S = dD and V = D, Eq. (26) becomes 



// 

J Jai 



XVX-hdS= I I I XV 2 XdV+ I I I \X7X\ 2 dV (27) 



IdD J J JD J J JD 

Since X (x) = for x e dD, 



II 

J Jai 



XV X -hdS = $ (2* 



'dD 

Also, from the PDE in (19), 



XV 2 XdV = -X / / / X 2 dV (29) 

D J J JD 

Substituting (28) and (29) into (27) gives 

° = - X J J J X 2 dV + J J J \VX\ 2 dV (30) 

For non-trivial solutions, X ^ at some points in D and hence by continuity of X, 
j j j D X 2 dV > 0. Thus (30) can be rearranged, 

fff n \VX\ 2 dV , x 

Since X is real, the the eigenvalue A is also real. 

If VX = for all points in D, then integrating and imposing the BC X (x) = 
for x G dD gives X = for all x e -D , i.e. the trivial solution. Thus VX is nonzero 
at some points in D, and hence by continuity of VX, j j j D \VX\ 2 dV > 0. Thus, 
from (31), A > 0. 
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6.3 Orthogonality of eigen-solutions to 
Sturm-Liouville problem 

Ref: Myint-U & Debnath §7.2 

Suppose vi, i> 2 are two eigenfunctions with eigenvalues Ai, A2 of the 3D Sturm- 
Liouville problem 

V 2 v + Xv = 0, x£D 
v = 0, x G dD 

We make use of Green's Formula (25) with V = D, S = dD, which holds for any 
functions v\, i> 2 defined on a volume V with smooth closed connected surface S. Since 
v 1 = = v 2 on dD and V 2 t>i = — Aii>i and V 2 f 2 = — A 2 t>2, then Eq. (25) becomes 




D 



= (Ai - A 2 ) / / / ^i^dV 

Thus if Ai ^ A 2 , 




Vl v 2 dV = 0, (32) 
and the eigenfunctions v 1 , v 2 are orthogonal. 



D 



7 Heat and Wave problems on a 2D rectangle, ho- 
mogeneous BCs 

Ref: Guenther & Lee §10.2, Myint-U & Debnath §9.4 - 9.6 
[Nov 7, 2006] 

7.1 Sturm-Liouville Problem on a 2D rectangle 

We now consider the special case where the subregion D is a rectangle 

D = {(x, y) : < x < x , < y < y } 
The Sturm-Liouville Problem (19) becomes 

^ + ^ + Xv = °» (x,y)£D, (33) 
v(0,y) = v(x ,y) = 0, < y < y , (34) 
v(x,0) = v(x,y ) = 0, < x < x . (35) 
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Note that in the PDE (33), A is positive a constant (we showed above that A had to 
be both constant and positive). We employ separation of variables again, this time 
in x and y: substituting v (x, y) = X (x) Y (y) into the PDE (33) and dividing by 
X (x) Y (y) gives 

Y" , X" 
Y +A= -X 

Since the l.h.s. depends only on y and the r.h.s. only depends on x, both sides must 
equal a constant, say /i, 

Y" X" 

^ + A = - x = /i (36) 
The BCs (34) and (35) imply 

X (0) Y (y) — X (x ) Y (y) — 0, < y < y , 

X (x) Y(0) = X (x) Y (y ) = 0, < x < x . 

To have a non-trivial solution, Y (y) must be nonzero for some y G [0,y ] and X (x) 
must be nonzero for some x G [0, x ], so that to satisfy the previous 2 equations, we 
must have 

X (0) = X (x ) = Y (0) = Y (y ) = (37) 
The problem for X (x) is the ID Sturm-Liouville problem 

X" + fiX = 0, < x < x (38) 
X(0) = X(x )=0 

We solved this problem in the chapter on the ID Heat Equation. We found that for 
non-trivial solutions, \i had to be positive and the solution is 

(TTX7[X \ ( TTX71 \ ^ 

J , IJ " m= \x~o~) ' "1=1,2,3,... (39) 

The problem for Y (y) is 

Y" + uY = 0, 0<y<y , (40) 
Y(0) = Y(y ) = 
where v — A — /i. The solutions are the same as those for (38), with v replacing ji: 

Y n (y) = b n sin(^^J , u n =(^j , n = 1,2,3,... (41) 
The eigen-solution of the 2D Sturm-Liouville problem (33) - (35) is 

(TYITCX \ ( TTTTy \ 
sin , m,n = 1,2,3, ... (42) 
xo J \ Uo / 

with eigenvalue 

/ 2 2 

2 / m n 



Amn l^m ~\~ ( 2 2 

■^0 ^0 
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7.2 Solution to heat equation on 2D rectangle 

The heat problem on the 2D rectangle is the special case of (7), 

u * = d^ + d^> t>0, 

u(x,y,t) = 0, (x,y) & dD, 
u(x,y,Q) = f(x,y), (x, y) G D, 

where D is the rectangle D = {(x,y) : < x < x , < y < y }- We reverse the 
separation of variables (9) and substitute solutions (14) and (42) to the T (t) problem 
(13) and the Sturm-Liouville problem (33) - (35), respectively, to obtain 

u mn (x,y,t) = A mn sm ( — — 1 sin I — J e w 

. .m'Kx\ . / rmy\ ~^(^+^\ 
= A mn sin sin e V x o y oJ 



%o J \ Do 

To satisfy the initial condition, we sum over all m, n to obtain the solution, in 
general form, 

oo oo oo oo ✓ , 

mux \ . I nixy \ _ x 



OO OO OO OO / 

u (x, y,t) = J2Yl Umn ( x ' y> *) = Yl Yl Amn sin ( ) sin ( — ) e~ w (43) 

m=ln=l ro=ln=l V ^0 / \ 2/0 / 

Setting £ = and imposing the initial condition -u (x, y, 0) = / (x, y) gives 

OOOO OOOO / \/ \ 

/ (x, y) =u (x, y, 0) = V] V] A mn w mn (x, J/) = Y] Y] sin ( ] sin ( — J 



where v mn (x,y) = sin ^^^j sin are the eigenfunctions of the 2D Sturm- 

Liouville problem on a rectangle, (33) - (35). Multiplying both sides by Vmh(x,y) 
(m, h — 1, 2, 3, ...) and integrating over the rectangle D gives 

/ / f(x,y) y) dA (44) 

J Jd m=ln=l J Jd 

where dA = dxdy. Note that 

f f f X ° ( 7TlTTx\ ( Thnx\ 

J J v mn (x, y) Vran (x, y) dA = J sin f —— J sin f — — J dx 

f yo ( nrry\ ( hiry\ 
x sin sin dy 

Jo V yo J V. yo J 

Xo2/o/4, m = m and n = h 
0, otherwise 
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Thus (44) becomes 

J J f( x , V) v mh (x, y) dA = ^j^x y 

Since to, h are dummy variables, we replace to by to and n by n, and rearrange to 
obtain 

A mn = / / f (x,y)Vrnn(x,y)dA 

xoVo J Jd 

4 r° f yo r . . (rrmx\ (nwy\ , , 

= / f(x, y) sin sin dydx (45) 

xoVo Jo Jo \ x J \ yo J 

7.3 Solution to wave equation on 2D rectangle, homogeneous 
BCs 

The solution to the wave equation on the 2D rectangle follows similarly. The general 
3D wave problem (8) becomes 

u tt = q^ 2 + W (x,y)eD, t>0, 

u(x,y,t) = 0, (x,y)edD, 
u(x,y,0) = f(x,y), (x, y) G D, 
u t (x,y,0) = g(x,y), (x, y) e D, 

where D is the rectangle D = {(x,y) : < x < x Q , < y < yo}- We reverse the 
separation of variables (9) and substitute solutions (16) and (42) to the T (t) problem 
(15) and the Sturm Liouville problem (33) - (35), respectively, to obtain 

u mn (x, y, t) = sin {^—^j sin cos (v^w) + Pnm sin (V^wt)) 

- sin - 



/ rmrx\ ( nny 
= sin sin 

V / V yo 



m 2 n 2 \ „ I m 2 n 2 



X I « nm COS I Trt\l — + ~2 + PnmSm Tit J — + 2 

x o Vo J \ \ x o Vo 

To satisfy the initial condition, we sum over all m, n to obtain the solution, in 
general form, 

(x, y,t) = ^^2 u mn (x, y, t) (46) 



oo oo 

U [ 

m=l n=l 

Setting t — and imposing the initial conditions 



u (x, y, 0) = / (x, y) , u t (x, y, 0) = g (x, y) 



11 



gives 

oo oo 

/ (x, y) = u (x, y, 0) = ^ ^2 a ™ v "™ ( x > v) 

m=l n=l 

v v -» / rmxx \ ( nny \ 
= > > «m« sin sin 

^ti V *o ; V yo J 

oo oo 

g (x, y) = u t (x, y,0) = ^2^2 \Z\^nl3 mn v m n (x, y) 

m= 

oo oo 

= ^ ^ V A mn /3; 



m=l n=l 



m7ra; \ / nny 
sin I sin 



where v mn (x,y) = sin sin (j^j are the eigenfunctions of the 2D Sturm Li- 

ouville problem on a rectangle, (33) - (35). As above, multiplying both sides by 
Vmn (x, y) (m, h = 1, 2, 3, ...) and integrating over the rectangle D gives 

a m n = / / / (x,y)v mn (x,y) dxdy 

x yo J J D 

Pmn = 7= / / g (x,y)v mn (x,y)dxdy 

XoyoV^mn J JD 

8 Heat and Wave equations on a 2D circle, homo- 
geneous BCs 

Ref: Guenther & Lee §10.2, Myint-U & Debnath §9.4, 9.13 (exercises) 

We now consider the special case where the subregion D is the unit circle (we may 
assume the circle has radius 1 by choosing the length scale I for the spatial coordinates 
as the original radius): 

D={(x,y):x 2 + y 2 <l} 
The Sturm-Liouville Problem (19) becomes 

d 2 v d 2 v 

x -2 + W +Xv = °> A ( 47 ) 

v(x,y) = 0, x 2 + y 2 = l, (48) 

where we already know A is positive and real. It is natural to introduce polar coor- 
dinates via the transformation 

x = rcos9, y = rsm6, w (r,9,t) = u (x,y,t) 

for 

< r < 1, -7i < 9 < 7i. 



12 



- d 2 v d 2 v 1 d ( dw\ 1 d 2 w 
V 2 r — ~ — — ~\~ — — — — -— r — + 



You can verify that 



dx 2 dy 2 r dr \ dr J r 2 d9 2 
The PDE becomes 

1 d ( dw\ 1 d 2 w , „ , 4 . 

;*( r *j + ^a9? +At " = ' 0<r<l,-,<9<, (49) 

The BC (48) requires 

w(l,6) = 0, -n<9<n. (50) 
We use separation of variables by substituting 

w (r, 6) = R (r) H (6) (51) 

into the PDE (49) and multiplying by r 2 / (R (r) if (0)) and then rearranging to obtain 

d ( dR\ 1 , 2 d 2 H 1 



dr \ dr J R (r) d6 2 H (9) 

Again, since the l.h.s. depends only on r and the r.h.s. on 9, both must be equal to 
a constant //, 

d / dR\ 1 , 9 d 2 H 1 



dr V r dr J R (r) + Ar " d^ (0) ~ (52) 
The BC (50) becomes 

w(l,9) = R(l)H(6) = 
which, in order to obtain non-trivial solutions (H {0) ^ for some 9), implies 

R(1) = (53) 

In the original (x, y) coordinates, it is assumed that v (x, y) is smooth (i.e. contin- 
uously differentiable) over the circle. When we change to polar coordinates, we need 
to introduce an extra condition to guarantee the smoothness of v (x,y), namely, that 

w (r, — 7r) = w (r, 7r) , wq (r, — ir) = we (r, ir) . (54) 

Substituting (51) gives 

H(-n)=H(n), ^ (-tt) = ^ (tt) . (55) 

The solution v (x, y) is also bounded on the circle, which implies R (r) must be 
bounded for < r < I. 
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The problem for H {6) is 

<I2H +HH(6) = Q ] H(-«) = H(*), ^ (-tt) = ^ (tt) . (56) 



d9 2 ^ y J ' v y v 7 ' d0 v y d# 

You can show that for /i < 0, we only get the trivial solution H (6) = 0. For (j, — 0, 
we have if (#) = const, which works. For /x > 0, non-trivial solutions are found only 
when jj, = m 2 , 

H m (#) = a m cos (m6) + b rn sin (m6) 

Thus, in general, we may assume A = m 2 , for m — 0, 1, 2, 3, ... 
The equation for i? (r) in (52) becomes 

d ( dR\ l,o o 
Rearranging gives 

r 2 ^ + r ^ + (\r 2 - m 2 ) R m = 0; R m (1) = 0, |fi ra (fl)|<oo (57) 
We know already that A > 0, so we can let 

s = V\r, R m (s) = R m (r) 

so that (57) becomes 



s 2 - 



2 + + (s 2 - m 2 ) R m = 0; R m (v^) = 0, |£ m (0)|<oo (58) 

The ODE is called Bessel's Equation which, for each m — 0, 1, 2, ... has two linearly 
independent solutions, J m (s) and y m (s), called the Bessel functions of the first and 
second kinds, respectively, of order m. The function J m (s) is bounded at s = 0; 
the function Y m (s) is unbounded at s = 0. The general solution to the ODE is 
R m (s) = c m iJ rn (s) + c m 2Y m (s) where c mn are constants of integration. Our bound- 
edness criterion \R m (0)| < oo at s = implies c 2m = 0. Thus 

R m (s) = c m J m (s) , R m (r) = c m J m (y/Xrj . 

[Nov 9, 2006] 

The Bessel Function J m (s) of the first kind of order m has power series 

(— l) fc s 2k + m 

Jm ^ = fci (A; + m)!22fc+™- (59) 



fc=0 



J m (s) can be expressed in many ways, see Handbook of Mathematical Functions by 
Abramowitz and Stegun, for tables, plots, and equations. The fact that J m (s) is 
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expressed as a power series is not a drawback. It is like sin (x) and cos (x), which are 
also associated with power series. Note that the power series (59) converges absolutely 
for all s > and converges uniformly on any closed set s G [0, L]. To see this, note 
that each term in the sum satisfies 



(-1)* 



k c 2k+m 



k\ (k + m)\2 2k+r ' 
Note that the sum of numbers 

£ 

k=0 



< 



j2k+r 



k\ (k + m)\2 2k + r 



^2k+r, 



k\ (k + m)\2 2k +^ 



converges by the Ratio Test, since the ratio of successive terms in the sum is 



£2(fc + l) + m 

(fc+l)!(fc+l+m)!2 2 ( fc + 1 )+'"' 



£2^ + 771 



k\{k+m)\2 2k + m 

Thus for k > N = \L/2] , 



< 



(fc + l)(fc + l+m)4 " (A: + l) 2 4 V 2 (^ + 1 ) 



£/ 2(fc+l) + n 



(fc+l)!(fc+l+m)!2 2 ( fc + 1 )+ r 



£2fc + r, 



< 



L 



2{N + l) 



< 1 



k\(k+m)\2 2k + m 

Since the upper bound is less than one and is independent of the summation index k, 
then by the Ratio test, the sum converges absolutely. By the Weirstrass M-Test, the 
infinite sum in (59) converges uniformly on [0, L\. Since L is arbitrary, the infinite 
sum in (59) converges uniformly on any closed subinterval [0, L] of the real axis. 

Each Bessel function J rn (s) has an infinite number of zeros (roots) for s > 0. Let 
j m ,n be the n'th zero of the function J m (s). Note that 

1 2 3 

J (s) .70,1 = 2.4048 j 0;2 = 5.5001 j 0f3 = 8.6537 
Ji(s) = 3.852 ji i2 = 7.016 j 1>3 = 10.173 

The second BC requires 

/? m (i) = R m (yX) = j m (yX) = o 

This has an infinite number of solutions, namely v^A = jm,n for n = 1, 2, 3, .... Thus 
the eigenvalues are 

A mn = 3m,ni rn,n = 1,2,3, ... 

with corresponding eigenfunctions J m {rj m , n ). The separable solutions are thus 



(x,y) 



(r,9) 



Jo(rj ,n) n = 1,2,3, ... 

J m (rj m ,n) cos (m6) + +(3 mn sin (m6>)) m, n = 1, 2, 3, ... 

(60) 
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8.1 Solution to heat equation on the 2D circle, homogeneous 
BCs 

The heat problem on the 2D circle is the special case of (7), 

(P'u d 2 u 

Ut = dx~^dy~ 2 ' (x '^ GjD ' t>0 ' 
u(x,y,t) = 0, (x,y)edD, 

u(x,y,Q) = f(x,y), (x, y) G D, 

where D is the circle D = {(x,y) : x 2 + y 2 < 1}. We reverse the separation of vari- 
ables (9) and substitute solutions (14) and (42) to the T (t) problem (13) and the 
Sturm Liouville problem (47) - (48), respectively, to obtain 

u mn (x, y, t) = v mn (x, y) e ~ Xmnt = v mn (x, y) e~ J ^ nt 

where v mn is given in (60). 

To satisfy the initial condition, we sum over all m, n to obtain the solution, in 
general form, 

oo oo oo oo 



u(x,y,t) = ^2^2u mn (x,y,t) = ^2^2v mn (x,y)e 

m=l n=l m=l n=l 

oo oo 

= ^5Zj m (rj min )(a mn cos(m^)+/3 mra sin(m^))e- w (61) 



m=l n=l 



where r = x 2 + y 2 and tan 6 = y/x. Setting t — and imposing the initial condition 

u (x, y, 0) = / (x, y) gives 

oo oo 

/ (x, y)=u (x, y, 0) = Yl Vmn ( x ' y) 

m=l n=l 

We can use orthogonality relations to find a mn and (5 mn . 

9 The Heat Problem on a square with inhomoge- 
neous BC 

We now consider the case of the heat problem on the 2D unit square 

D = {(x,y):0<x,y<l}, (62) 
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where a hot spot exists on the left side, 



u t = V u, (x, y) G D 

( i \ ( u /e, {x = 0, \y-y \ <e/2} 
u(x,y,t) = < n (63) 
[ 0, otherwise on c>L> 

u(x,j/,0) = 

where the hot spot is confined to the left side: < yo — e/2 < y < y + e/2 < 1. 

9.1 Equilibrium solution and Laplace's Eq. on a rectangle 

Ref: Guenther & Lee §8.1, 8.2, Myint-U & Debnath §6.6, 8.7 

As in the ID case, we first find the equilibrium solution ue (x,y), which satisfies 
the PDE and the BCs, 



V 2 u E = 0, (x,y)eD 

u /e {x = 0, |y - y \ < e/2} 
otherwise on dD 



ue (x,y) 



The PDE for ue is called Laplace's equation. Laplace's equation is an example of 
an elliptic PDE. The wave equation is an example of a hyperbolic PDE. The heat 
equation is a parabolic PDE. These are the three types of second order (i.e. involving 
double derivatives) PDEs: elliptic, hyperbolic and parabolic. 

We proceed via separation of variables: ue (x,y) = X (x) Y (y), so that the PDE 
becomes ^ y// 

where A is constant since the l.h.s. depends only on x and the middle only on y. The 
BCs are 

y(0) = Y(l) = 0, X(1) = 

and 

\ otherwise 
We first solve for Y (y), since we have 2 easy BCs: 

Y" + AF = 0; F(0) = F(1)=0 

The non-trivial solutions, as we have found before, are Y n = sin (mry) with \ n = n 2 ir 2 , 
for each n = 1, 2, 3, ... Now we consider X (x): 

X" - nVX = 
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and hence 

X (x) = Cl e nwx + c 2 e- n ™ 
An equivalent and more convenient way to write this is 

X (x) = C3 sinhn7r (1 — x) + C4 cosher (1 — x) 

Imposing the BC at x — 1 gives 

X(l) = c 4 = 

and hence 

X (x) = C3 sinh nir (1 — x) 

Thus the equilibrium solution to this point is 

00 

ue (x, y) = A n sinh (nir (1 — x)) sin (niry) 

n=l 

You can check that this satisfies the BCs on x = 1 and y — 0, 1. Also, from the BC 
on x = 0, we have 

Mb (0, y) = Y] A n sinh (rar) sin (ran/) = J ' y 2/0 ' < E ^ 

^— f otherwise 

71=1 ^ 

Multiplying both sides by sin (miry) an integrating in y gives 

» /•! rvo+e/2 Uq 

y A n sinh (rar) / sin (ran/) sin (miry) dy — — sin (miry) dy 

n= l JO Jyo-e/2 £ 



From the orthogonality of sin's, we have 

1 

2 Jyo-e/2 £ 



1 rvo+e/2 

A m sinh (rrar) - = / — sin (miry) dy 



Thus, 



2u ry»+ £ / 2 
A ™ = — r-TT r / sin (miry) dy 



VQ+e/2 



esinh (mir) J yo - £ / 2 



2uo cos (miry) 

e sinh (mir) 
2u 



mir 



emir sinh (mir 
4-u sin (miry ) sin ( !Z y i ) 
em7r sinh (mir) 
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J/o-e/2 

cos (rrar (y — e/2)) — cos (mir (y + s/2))) 



Thus 



( , 4M ^ sin (n7ry )sin(^f ) . 

We (x, v) = > — ; Sinn l™ 71 " (1 — x)) sin (ran/) 

y ,yj en ^ nsinh wr v v )> \ y) 

n=l v ' 



To solve the transient problem, we proceed as in 1-D by defining the function 

v (x, y,t) = u (x, y, t) - u E (x, y) 

so that v (x, y, t) satisfies 

v t = V 2 v 
v = on dD 
v(x,y,0) = f (x,y) -u E (x,y) 

9.2 First term approximation 

To approximate the equilibrium solution u E (x,y), note that 

sinhnvr (1 - x) _ e n < x ~^ - e'^ 1 "^ 
sinh mr 

For sufficiently large n, we have 

sinhn7r (f — x) e ™r(l-s) 



sinh nn 



Thus the terms decrease in magnitude (x > 0) and hence u E (x,y) can be approxi- 
mated the first term in the series, 

4g, S in(^)Bin(f) 
v ' y; £7r sinh(7r) v v v y> 

A plot of sinh (jr (1 — x)) sin (V?/) is given below. The temperature in the center of 
the square is approximately 

/I 1\ 4^ sm(^)sm(f) 
V 2 2/ sinh(7r) V2/ V2 



9.3 Easy way to find steady-state temperature at center 

For y Q = 1/2 and e — 1, we have 

1 1\ 4w sinh(^) u 



2' 2/ vt sinh(7r) 4 
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X 



Figure 1: Plot of sinh (ir(l — x)) sm(ny). 



It turns out there is a much easier way to derive this last result. Consider a plate 
with BCs u = uo on one side, and u = on the other 3 sides. Let a = ue (|, |). 
Rotating the plate by 90° will not alter u E (|, ^), since this is the center of the plate. 
Let tiEsum be the sums of the solutions corresponding to the BC u = u on each of the 
four different sides. Then by linearity, iiEsum = «o on all sides and hence UEsum = «o 
across the plate. Thus 

MO = UEsum Q, ^ = 4« 

Hence a = ue (|, |) = «o/4. 



9.4 Placement of hot spot for hottest steady-state center 

Note that 

4w sin (7n/ ) sin (f) 



u E 



1 1 

2' 2 



sinh 



(D 



£7r sinh (it) 
4m sinh (|) sin (iry ) sin ( 
7r sinh f 7r) 



— sin (7ry ) - 

4 7T 
MO . / v 2 

— sin (ny ) - 

4 7T 



e 
(f ) 

2 



2tt 



sin (7tj/ ) 
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for small e. Thus the steady-state center temperature is hottest when the hot spot is 
placed in the center of the side, i.e. y = 1/2. 



9.5 Solution to inhomogeneous heat problem on square 

We now use the standard trick and solve the inhomogeneous Heat Problem (63), 

u t = V 2 m, (x,y) e D 

»<*,»«) = \ u f ^=^r*i<f ,64) 

I 0, otherwise on oD 

u(x,y,0) = f(x,y) 

using the equilibrium solution u E . We define the transient part of the solution as 

v (x, y,t) = u (x, y, t) - u E (x, y) 

where u (x,y,t) is the solution to (64). The problem for v (x,y,t) is therefore 

v t = V 2 v, (x,y) G D 
v(x,y,t) = 0, (x,y)edD 
v(x,y,0) = f(x,y)-u E (x,y) 

This is the Heat Problem with homogeneous PDE and BCs. We found the solution 
to this problem above in (43), (provided we use x = 1 = y in (43); don't mix up 
the side length y in (43) with the location y of the hot spot in the homogeneous 
problem): 

oo oo 

v (x, y,t) = ^ A mn sin (mirx) sin (mcy) e - 7r2 ( m2 +" 2 )* 

m=l n=l 

where 

1 rl 

(/ (x, y) — ue (x, y)) sin (mux) sin (niry) dydx 



'o J o 

Thus we have found the full solution u (x, y, t). 



10 Heat problem on a circle with inhomogeneous 

BC 

Consider the heat problem 

u t = V 2 m, (x,y) e D 

u(x,y,t) = g(x,y), (x,y)edD (65) 

u(x,y,0) = f(x,y) 

21 



where D = {(x, y) : x 2 + y 2 < 1} is a circle (disc) of radius 1. 



10.1 Equilibrium solution and Laplace's Eq. on a circle 

Ref: Guenther & Lee §8.1, 8.2, Myint-U & Debnath §8.4, 8.5 

To solve the problem, we must first introduce the steady-state u = ue (x,y) which 
satisfies the PDE and BCs, 

V 2 u E = 0, (x,y)eD, 
u E (x,y) = g(x,y), (x,y)edD. 

As before, switch to polar coordinates via 

x = rcos9, y = rsin9, we (r,9) — ue (x,y) 

for 

< r < 1, -7i < 9 < 7T . 
The problem for u E becomes 

1 d ( dw E \ 1 d 2 w E „ , s 

r*( r ^) + ^ = °- <><r<l, -»<«<, (66) 

w (r, — 7r) = u> (r, 7r) , we (r, — tt) = we (r, tt) , (67) 

|iu (0, 0) | < oo (68) 

w(l,9) = g(9) , -n<9<n (69) 

where g (9) = g (x,y) for (x,y) G <9D and 9 = arctan (y/x). 
We separate variables 

w ( r , 9) = R (r) if (0) 

and the PDE becomes 

r d ( dR\ d 2 H 1 



R (r) rfr V / dB 2 H (9) 

Since the l.h.s. depends only on r and the r.h.s. on 9, both must be equal to a 
constant /i, 

r d / diA d 2 if 1 

r ~n = - J7v>. tt 7ns = V 



R (r) dr V dr y if (9) 

The problem for if (9) is, as before, 

d 2 ff tt , n \ tt / \ tt i \ dH , . dH , . 
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with eigen-solutions 

H m (9) = a m cos (m6) + b m sin (m9) , to = 0, 1, 2, ... 

The problem for R (r) is 

d ( dR\ ~ 2 d 2 R dR 2 

= r— r— — m R = r —— + r— to R 

dr \ dr J dr A dr 

Try R (r) = r a to obtain the auxiliary equation 

a (a — 1) + a — m 2 = 

whose solutions are a = ±m. Thus for each m, the solution is R m (r) = c\r m + C2r~ m . 
For m > 0, r~ m blows up as r — > 0. Our boundedness criterion (68) implies c 2 = 0. 
Hence 

(r) = c m r m 

where c m are constants to be found by imposing the BC (69). The separable solutions 
satisfying the PDE (66) and conditions (67) to (68), are 

w m (r, 9) = r m (Am cos (m6) + B m sin (mO)) , m — 0, 1, 2, ... 

The full solution is the infinite sum of these over to, 

oo 

u E (x, y) = w E (r, 9) = ^ r m (A m cos (md) + B m sin (mO)) (70) 

m=0 

We still need to find the A m , B m . 
Imposing the BC (69) gives 

oo 

(A m cos (m9) + B m sin (m9)) = g (6) (71) 

m=0 

The orthogonality relations are 



/' 



cos (m9) sin (nO) d6 = 



M cosM)cos(^) 1^ = f tt, TO = n (m>Q) ^ 
,r I sin (to6*) sin (n^) j j 0, to 7^ n 

Multiplying (71) by sinn^ or cosnO and applying these orthogonality relations gives 

A = 1-J* g(6)d6 

A m = - f g(6)cos(m6)d6 (72) 
B m = - I g(9)sm(m9)d9 

7T i-TT 
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Figure 2: Setup for hot spot problem on circle. 



10.1.1 Hot spot on boundary 

Suppose 

HO) = 



A -7T < 9 < 9 



otherwise 

which models a hot spot on the boundary. The Fourier coefficients in (72) are thus 

M = — , A m = — — — rsm mfc) 

27r mn [Vo + 7TJ 

B m = — — — - (cos (m6 ) - (-l) m ) 

mir (Vo + 7TJ 

Thus the steady-state solution is 

u ° , V 1 m / «o sin (m# ) , m w (cos (m# ) - (~l) m ) . , _A 
we = h > r — cos m& — n : sin (mO) 

m=l \ \ / \ / / 



10.1.2 Interpretation 

The convergence of the infinite series is rapid if r <C 1. If r ~ 1, many terms are 
required for accuracy. 

The center temperature (r = 0) at equilibrium (steady-state) is 

u E (0, 0) = w E (0,6) = ^ = i f £ (5) d# 

27T 27T / 
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i.e., the mean temperature of the circumference. This is a special case of the Mean 
Value Property of solutions to Laplace's Equation V 2 w = 0. 

We now consider plots of ue (x,y) for some interesting cases. We draw the level 
curves (isotherms) ue = const as solid lines. Recall from vector calculus that the 
gradient of %, denoted by Vue, is perpendicular to the level curves. Recall also 
from the physics that the flux of heat is proportional to Vu E - Thus heat flows along 
the lines parallel to Vue- Note that the heat flows even though the temperature is in 
steady-state. It is just that the temperature itself at any given point does not change. 
We call these lines the "heat flow lines" or the "orthogonal trajectories", and draw 
these as dashed lines in the figure below. 

Note that lines of symmetry correspond to (heat) flow lines. To see this, let n; 
be the normal to a line of symmetry. Then the flux at a point on the line is Vw ■ n^. 
Rotate the image about the line of symmetry. The arrow for the normal to the line of 
symmetry is now pointing in the opposite direction, i.e. — n/, and the flux is — Vw-nj. 
But since the solution is the same, the flux across the line must still be Vw • n/. Thus 



which implies Vw-nz = 0. Thus there is no flux across lines of symmetry. Equivalently, 
Vw is perpendicular to the normal to the lines of symmetry, and hence Vw is parallel 
to the lines of symmetry. Thus the lines of symmetry are flow lines. Identifying the 
lines of symmetry help draw the level curves, which are perpendicular to the flow 
lines. Also, lines of symmetry can be thought of as an insulating boundary, since 
Vm • n ( = 0. 



Use the BCs for the boundary. Note that the solution is symmetric with respect 
to the y-axis (i.e. even in x). The solution is discontinuous at {y = 0,x = ±1}, or 
{r = 1, 6 — 0, 7r}. See plot. 

(ii) —7i < 9q < — tt/2. The sum for ue is messy, so we use intuition. We start 
with the boundary conditions and use continuity in the interior of the plate to obtain 
a qualitative idea of the level curves and heat flow lines. See plot. 

(iii) 6q — > — 7r + (a heat spot). Again, use intuition to obtain a qualitative sketch 
of the level curves and heat flow lines. Note that the temperature at the hot point 
is infinite. See plot. The heat flux lines must go from a point of hot temperature to 
a point of low temperature. Since the temperature along boundary is zero except at 



Vw • ni = — Vw ■ 



(i) #o = 0. Then 
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e n < -7t/2= 



Figure 3: Plots of steady-state temperature due to a hot segment on boundary. 

9 = — 7r, then all the heat flux lines must leave the hot spot and extend to various 
points on the boundary 

10.2 Solution to inhomogeneous heat problem on circle 

We now use the standard trick and solve the inhomogeneous Heat Problem (65), 

u t = V 2 u, (x,y) e D 
u(x,y,t) = g(x,y), (x,y)edD (73) 
u(x,y,0) = f(x,y) 

using the equilibrium solution ue- We define the transient part of the solution as 

v (x, y } i) = u (x, y, t) - u E (x, y) 

where u (x,y,t) is the solution to (73). The problem for v (x,y,t) is therefore 

v t = V 2 f, (x,y) E D 
v(?,y,t) = 0, (x,y)edD 
v(x,y,0) = f (x,y) - u E (x,y) 

This is the Heat Problem with homogeneous PDE and BCs. We found the solution 
to this problem above in (61), 

oo oo 

v (x, y, t) = ^2 ^2 Jm ( r 3m,n) («mn cos (m9) + P mn sin (mO)) e ' Xmnt 

m=l n=l 

where a mn and (3 mn are found from orthogonality relations and f (x,y) — u E (x,y). 
Thus we have found the full solution u (x, y, t). 
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11 Mean Value Property for Laplace's Eq. 

Ref: Guenther & Lee §8.4, Myint-U & Debnath §8.4 

Theorem [Mean Value Property] Suppose v (x, y) satisfies Laplace's equation 
in a 2D domain D, 

V 2 v = 0, (x,y)ED. (74) 

Then at any point (x ,y ) in D, v equals the mean value of the temperature around 
any circle centered at (xo,yo) and contained in D, 

i r 

v(x ,y ) = — / v (xq + Rcos9,y + Rsin9) dO. (75) 

Note that the curve {(x + Rcos9,y + Rsin9) : — ir < 9 < n} traces the circle of 
radius R centered at (xo,yo)- 

Proof: To prove the Mean Value Property, we first consider Laplace's equation 
(74) on the unit circle centered at the origin (x,y) = (0,0). We already solved this 
problem, above, when we solved for the steady-state temperature ue that took the 
value g(9) on the boundary. The solution is Eqs. (70) and (72). Setting r = in 
(70) gives the center value 

u E (0,0)=A = ^-J g(9)d9 (76) 

On the boundary of the circle (of radius 1), (x,y) = (cos 9, sin 9) and the BC implies 
that u E = g (9) on that boundary. Thus, g (9) = u E (cos 9, sin 9) and (76) becomes 

i r 

u E (0,0) = A = — u E (cos 9, sin 9) d9. (77) 
We now consider the region 

B (xo, yo ) ( R ) = {( x ,y) : ( x - x o) 2 + (y- yof <R 2 }, 

which is a disc of radius R centered at (x,y) = (xo,yo)- Since B^ xo>yo ^(R) C D, 
Laplace's equation (74) holds on this disc. Thus 

V 2 v = 0, (x,y)eB {xo , yo) (R). (78) 
We make the change of variable 

X= X j^ , y = V R V() > u E {x,y) =v(x,y) (79) 

to map the circle B^ XO;yo ^ (R) into the unit disc {(x, y) : x 2 + y 2 < 1}. Laplace's equa- 
tion (78) becomes 

VVe = 0, (x,y) E {(x,y) : x 2 + y 2 < l} 
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where V 2 = (d 2 / dx 2 , d 2 / dy 2 ) . The solution is given by Eqs. (70) and (72), and we 
found the center value above in Eq. (77). Reversing the change of variable (79) in 
Eq. (77) gives 

1 r 

v(x ,y ) = — / v(x + Rcos9,y + Rsm9)d9 

as required. ■ 

For the heat equation, the Mean Value Property implies the equilibrium tempera- 
ture at any point (xq, yo) in D equals the mean value of the temperature around any 
disc centered at (xo,yo) and contained in D. 

12 Maximum Principle for Laplace's Eq. 

Ref: Guenther & Lee §8.4, Myint-U & Debnath §8.2 

Theorem [Maximum Principle] Suppose v (x, y) satisfies Laplace's equation 
in a 2D domain D, 

V 2 v = 0, (x,y)eD. 

Then the function v takes its maximum and minimum on the boundary of D, 3D. 

Proof: Let {xo,yo) be an interior point in D where v takes it's maximum. In 
particular, (x,y) is not on the boundary 3D of D. The Mean Value Property implies 
that for any disc B^ Xom ^{R) of radius R > centered at (xo,yo), 

1 r 

v(xo,yo) = — / v (xo + Rcos9,y + RsinO) d6 
= Average of v on boundary of circle 

But the average is always between the minimum and maximum. Thus v (xq, yo) must 
be between the minimum and maximum value of v (x, y) on the boundary of the disc. 
But since v (xq, yo) is the maximum of v on D, then the entire boundary of the circle 
must have value v (xq, yo). Since R > is arbitrary, this holds for all values in the disc 
B(x ,y ){R)- Keep increasing R until the disc hits the boundary of D. Then v (x ,yo) 
is a value of v along the boundary of D. 

A similar argument holds for the minimum of v. ■ 

Corollary If v = everywhere on the boundary, then v must be zero at every- 
where in D. 

For the homogeneous heat equation (i.e. no sinks/sources), this implies the equi- 
librium temperature attains its maximum/minimum on the boundary. 
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If the boundary is completely insulated, Vi> = on dD, then the equilibrium 
temperature is constant. To see this, we apply result (26) to the case where v is the 
solution to Laplace's equation with Vi> = on dD, 

0= I [ vVv-hdS= [[[ (vV 2 v + \Vv\ 2 ) dV = [ [ [ \Vv\ 2 dV 
J JdD J J Jd J J Jd 

Thus, by continuity, |Vi> | = and Vi> = everywhere in D and hence v = const in 

D. 



13 Eigenvalues on different domains 

In this section we revisit the Sturm-Liouville problem 

V 2 v + Xv = 0, x£D 
v = 0, x G dD 

and consider the effect of the shape of the domain on the eigenvalues. 
Definition Rayleigh Quotient 

, , f f fn • VvdV , 

Theorem Given a domain D C R 3 and any function v that is piecewise smooth 
on D, non-zero at some points on the interior of D, and zero on all of dD, then the 
smallest eigenvalue of the Laplacian on D satisfies 

A < R(v) 

and R(h) = A if and only if h (x) is an eigen-solution of the Sturm Liouville problem 
on D. 

Sketch Proof: We use result (26) derived for any smooth function v defined on 
a volume V with closed smooth surface S. 

If vVv-hdS= I I I (vV 2 v + Vv-Vv)dV 
J JdD J J Jd 

In the statement of the theorem, we assumed that v = on dD, and hence 

Vv ■ VvdV = - [ [ [ vV 2 vdV (81) 
d J J Jd 

Let {4> n } be an orthonormal basis of eigenfunctions on D, i.e. all the functions <p n 

which satisfy 

V 2 0n + A n 0„ = 0, xe£> 
4> n = 0, x G dD 
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and 

/ / / <f>n<f>mdV = \ lj m 71 

J J Jd { 0, n 

We can expand v in the eigenfunctions, 

oo 

V (x) = ^ A n$n (x) 
n=l 

where the A n are constants. Assuming we can differentiate the sum termwise, we 
have 

oo oo 

V 2 V = Ay 2 <Pn = -J2 X nAn<Pn (82) 

rt=l n=l 

The orthonormality property (i.e., the orthogonality property with f f f D (findV = 1) 
implies 

oooo ««« oo 

/ / / v 2 dV = J2J2 A " A ™ / / <Pn<PmdV = J2 A l (83) 

J J JD n =lm=l J J JD n=1 

and, from (82), 

K 2 vdV = - jr J2 \ n A n A m I!! ^ m dV = -J2 \ n A\ (84) 

m=l n=l J J JD n=1 

Substituting (84) into (81) gives 

/ / / VV ' VvdV = ~ 1 1 1 vV2ydV = £ XnA ™ ^ 
Substituting (83) and (85) into (80) gives 

nf D vvvvdv _ j:z 1 \nAi 




R v) = 



We assume the eigenfunctions are arranged in increasing order. In particular, A n > Ai. 
Thus 

V°° A, A 2 V°° A 2 

n{V) > V oo , 2 - A l V oo A2 ~ A l- 
Z-m=l n Z^n=l 

If v is an eigenfunction with eigenvalue Ai, then 



R(v) = ^ = X l 



If v is not an eigenfunction corresponding to Ai, then there exists an n > 1 such that 
A n > Ai and A n ^ 0, so that R (v) > X±. ■ 
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Theorem If two domains D and D in M 2 satisfy 

D g D i.e., D C D but D ^ D, 

then the smallest eigenvalues of the Sturm-Liouville problems on D and D, X 1 and 
Ai, respectively, satisfy 

Ai < Ai 

In other words, the domain D that contains the sub-domain D is associated with a 
smaller eigenvalue. 

Proof: Note that the Sturm-Liouville problems are 

V 2 v + Aw = 0, (x,y) e D 
v = 0, {x,y)edD 



V 2 v + Xv = 0, (x,y) G D 
v = 0, G 95 

Let i>i be the eigenfunction corresponding to Ai on D. Then, as we have proven 
before, 

\i = R{v x ), (86) 

where R(v) is the Rayleigh Quotient. Extend the function v\ continuously from D 
to D to obtain a function v\ on D which satisfies 




v 1 (x,y), (x,y)ED 

(x,y)eD, (x,y)<£D 



The extension is continuous, since V\ is zero on the boundary of D. Applying the 
previous theorem to the region D and function v\ (which satisfies all the requirements 
of the theorem) gives 

\i<R(vi). 

Equality happens only if v 1 is the eigenfunction corresponding to Ai. 

Useful fact [stated without proof]: the eigenfunction(s) corresponding to the 
smallest eigenvalue A\ on D are nonzero in the interior of D. 

From the useful fact, v\ cannot be an eigenfunction corresponding to Ai on D, 
since it is zero in the interior of D (outside D). Thus, as the previous theorem states, 

Ai <R (i)i). (87) 
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Since v 1 = outside D, the integrals over D in the Rayleigh quotient reduce to 
integrals over D, where Vi = v±, and hence 

R(v 1 )=R(v 1 ). (88) 

Combining (86), (87), and (88) gives the result, 

Ai < Ai. 



Example: Consider two regions, D\ is a rectangle of length xo, height yo and D 2 
is a circle of radius R. Recall that the smallest eigenvalue on the rectangle D x is 

2 ( 1 1 

An = 7T I -o + ~o 

The smallest eigenvalue on the circle of radius 1 is Aoi = Jqi w here the first zero 
of the Bessel function J (s) of the first kind is Jo,i = 2.4048. Since v^A multiplied 
r in the Bessel function, then for a circle of radius R, we'd rescale by the change of 
variable f = r/R, so that J m (^/Xr^j = J m (^^J~-^r^ where f goes from to 1. Thus 
on the circle of radius R, the smallest eigenvalue is 

Aoi = (^j , Jo,i = 2.4048. 

Suppose the rectangle is actually a square of side length 2R. Then 

7T 2 _ 4.934 fJo,iY 5.7831 

All — 777^ — — F^ - ' Aoi 



2R 2 R 2 ' Ui V R J R 2 

Thus, An < A i, which confirms the second theorem, since D 2 C D 1 , i.e., the circle 
is contained inside the square. 
Now consider the function 

You can show that 

Vv ■ Vv = 



dv x z 



dr 
and 

(-R (dv\ 2 



Jf D Vv.VvdA _ f:j«(%yrdrd9 _ 6 
IL 2 ^ j:j R v 2 rdrd9 R 2 



This confirms the first theorem, since f (r) is smooth on D2, « (-R) = (zero on the 
boundary of D2), and v is nonzero in the interior. 
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13.1 Faber-Kahn inequality 

Thinking about the heat problem on a 2D plate, what shape of plate will cool the 
slowest? It is a geometrical fact that of all shapes of equal area, the circle (disc) has 
the smallest circumference. Thus, on physical grounds, we expect the circle to cool 
the slowest. Faber and Kahn proved this in the 1920s. 

Faber-Kahn inequality For all domains D C M 2 of equal area, the disc has the 
smallest first eigenvalue Ai. 

Example. Consider the circle of radius 1 and the square of side length y/n. Then 
both the square and circle have the same area. The first eigenvalues for the square 
and circle are, respectively, 

Aisq = 7T 2 (- + -)= 2tt = 6.28, X WIRC = (J ,i) 2 = 5.7831 

\7l 71 J 

and hence \\sq > ^icirc, as the Faber-Kahn inequality states. 

14 Nodal lines 

Consider the Sturm-Liouville problem 

X7 2 v + \v = 0, xeD (89) 
v = 0, x G 3D 

Nodal lines are the curves where the eigenfunctions of the Sturm-Liouville problem 
are zero. For the solution to the vibrating membrane problem, the normal modes 
t) are zero on the nodal lines, for all time. These are like nodes on the ID 
string. Here we consider the nodal lines for the square and the disc (circle). 



14.1 Nodal lines for the square 

For the square, the eigenfunctions and eigenvalues are a special case of those we found 
for the rectangle, with side length xo = yo = a, 

/mrcx\ /mny\ n 2 , 2 2 , 

v mn {x,y) = sin ^— — J sin \— — J , \ mn = — (m + n ) , n, m = 1, 2, 3, ... 

The nodal lines are the lines on which v mn (x, y) = 0, and are 
ka la 

x = — , y — — , 1 < k < m — 1, 1 < I < n — 1 
m n 

for to, n > 2. Note that Vu (x,y) has no nodal lines on the interior - it is only zero 
on boundary dD. Since \ mn = \ nm , then the function f nm = Av rnn + Bv nm is also an 
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eigenfunction with eigenvalue A mn , for any constants A, B. The nodal lines for f nm 
can be quite interesting. 

Examples: we draw the nodal lines on the interior and also the lines around the 
boundary, where v nm = 0. 
(i) m — 1, n — 1. 



r ; i = ii ( — j sin ^ — j 



f 71X * 

V a 

This is positive on the interior and zero on the boundary. Thus the nodal lines are 
simply the square boundary dD. 
(ii) m — 1, n — 2. 

'ttx\ . / 2iry s 



V\2 = sin ^ — j sin 

The nodal lines are the boundary <9.D and the horizontal line y = a/2. 

(iii) m — 1, n — 3 

/7TX\ / 37Ty 

Di3 = sin — sin 

\ o / \ a 

The nodal lines are the boundary dD and the horizontal lines y = a/3, 2a/ 3. 

(iv) to = 3, n — 1 

V3i = sm {—) 8111 ItJ 

The nodal lines are the boundary dD and the vertical lines x = a/3, 2a/3. 

(v) Consider t>i 3 — v 3 i. Since A 3 i = Ai 3 = 107r 2 /a 2 , this is a solution to SL 
problem (89) on D with A = A i3 . To find the nodal lines, we use the identity 
sin 30 = (sin 9) (3 - 4 sin 2 9) to write 

„ 13 = sin (^) sin (^) (3-4sin 2 (^)) 

/ttx\ . /?ry\ / . 2 /7TX\ . 2 /ny 

I I sin I I I sin I 1 — sin I 

a 



V13 — V31 = 4 sin ^ — J sin ^ — J |^sin ^ — J — sin 



The nodal lines are the boundary of the square, dD, and lines such that 



I I — sin I — I = I sin I I — sin I — I I I sin I I -I- sin I — I I 



= siir ' tJ - sin ItJ = l sin ItJ ~ sin ItJ ) l sin ItJ + sin v „ 



i.e., 

sin I ( 

a / V a ) 



. nx\ . ( , ny\ 
sin — = sin ± — 



Thus 



— = ± h K7T, any integer k. 

a a 
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Hence the nodal lines are 

y = ±x + ka 

for all integers k. We're only concerned with the nodal lines that intersect the interior 
of the square plate: 

y = x, y = —x + a 

Thus the nodal lines of v 13 — v 31 are the sides and diagonals of the square. 

Let Dt be the isosceles right triangle whose hypotenuse lies on the bottom hor- 
izontal side of the square. The function t> 13 — v 31 is zero on the boundary dD?, 
positive on the interior of D?, and is thus the eigenfunction corresponding to the first 
(smallest) eigenvalue A = A13 of the Sturm-Liouville problem (89) on the triangle Dt- 
In this case, we found the eigenfunction without using separation of variables, which 
would have been complicated on the triangle. 

(vi) With v 13, i>3i given above, adding gives 

v 13 + „ 31 = 4sin (=) S i„ (!?){»_ ( sin * + sm * (^)) } 

The nodal lines for t>i 3 + t> 3 i are thus the square boundary and the closed nodal line 
defined by 




Let D c be the region contained within this closed nodal line. The function — (v 13 + V31) 
is zero on the boundary dD c , positive on the interior of D c , and thus — (v 13 + V31) is 
the eigenfunction corresponding to the first (smallest) eigenvalue A i3 of the Sturm- 
Liouville problem (89) on D c . 

(vii) Find the first eigenvalue on the right triangle 

{o<y<V3x, 0<x<l}. 

Note that separation of variables is ugly, because you'd have to impose the BC 

X (x) Y (VSx^j = 

We proceed by placing the triangle inside a rectangle of horizontal and vertical side 
lengths 1 and v^3 ? respectively. The sides of the rectangle coincide with the perpen- 
dicular sides of the triangle. Thus, all eigenfunctions v mn for the rectangle are already 
zero on two sides of the triangle. However, any particular eigenfunction v mn will not 
be zero on the triangle's hypotenuse, since all the nodal lines of v mn are horizontal 
or vertical. Thus we need to add multiple eigenfunctions. However, to satisfy the 
Sturm-Liouville problem, all the eigenfunctions must be associated with the same 
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eigenvalue. For the rectangle with side lengths and 1, the eigenvalues are given 
by 




We create a table of 3m 2 + n 2 and look for repeated values: 



n, m 


1 


2 


3 


4 


1 


4 


13 


28 




2 


7 


16 


31 




3 


12 


21 


36 




4 


19 


28 


43 




5 


28 


37 







The smallest repeated value of 3m 2 + n 2 is 28: 

A31 = A 2 4 = A15 = 287r 2 /3 

The linear combination of the corresponding eigenfunctions is itself an eigenfunction 
of the SL problem (89) on the rectangle D, 

V = Av 31 + Bv 2 4 + CV15. 

We wish v to also be zero along the hypotenuse of D T 2, and seek A, B, C which 
satisfy 

AV 31 + Bv 24 + Cv 15 = 

on y = V3x. If we can find such constants A, B, C then we have found an eigenfunc- 
tion of the SL problem (89) on the triangle Dt2- If this eigenfunction is positive on 
the interior of D?2 then we have also found the smallest eigenvalue (in this case A31). 

If we cannot find such A, B, C, then we seek the next largest repeated value 
of 3m 2 + n 2 and continue as above. If we can find (A, B, C) such that v = on 
y = \^3x, but v is zero on the interior of Dt2, then A31 is not the smallest eigenvalue 
and we must pursue another method of finding the smallest eigenvalue. The Rayleigh 
quotient can help us identify an upper bound. 

14.2 Nodal lines for the disc (circle) 

For the disc of radius 1, we found the eigenfunctions and eigenvalues to be 

v m ns = Jm (rj m ,n) smm9, v mnC = J m (rj m>n ) cos m6> 
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with 

Amn = ^ 2 3m,ni 71,771 = 1,2,3,... 

The nodal lines are the lines on which v mn c = or v mn s = 0. 
Examples. 

(i) m = 0, n = 1. 

voi = Jo (rjo,i) 

The nodal lines are the boundary of the disc (circle of radius 1). 

(ii) m = 0, n = 2. 

V02 = Jo (rj , 2 ) 

The nodal lines are two concentric circles, one of radius r = 1, the other or radius 

r = Jo,i/Jo,2 < 1- 

(hi) m = 1 , n — 1 and sine. 

vus = Jo(rji,i)sm9 

The nodal lines are the boundary (circle of radius 1) and the line 9 = — 7r, 0, ir 
(horizontal diameter). 



15 Steady-state temperature in a 3D cylinder 

Suppose a 3D cylinder of radius a and height L has temperature u (r,9, z,t). We 
assume the axis of the cylinder is on the z-axis and (r, 9, z) are cylindrical coordinates. 
Initially, the temperature is u (r, 9, z, 0). The ends are kept at a temperature of u — 
and sides kept at u (a, 9, z,t) = g (9, z). The steady-state temperature ue (r, 9, z) in 
the 3D cylinder is given by 

VVe = 0, -7i < 9 < 7i, < r < a, < z < L, (90) 

u E (r,9,0) = u E (r,9,L) = 0, -tt < 9 < tt, < r < a, (91) 
u E (a,9,z) = g(9,z), -n < 9 < n, < z < L. 

In cylindrical coordinates (r, 0, z), the Laplacian operator becomes 

„, I d f dv\ 1 (9 2 w 9 2 v 

V v = r — H 1 

r dr \ dr J r 2 d9 2 dz 2 

We separate variables as 

v(r,9,z) = R(r)H(0)Z(z) 
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so that (90) becomes 



1 d ( dR(r)\ 1 d 2 H{0) 1 d 2 Z{z) _ 



rR (r) rfr V / r 2 #(#) d6 2 Z (z) dz 2 

Rearranging gives 

1 d ( dR(r)\ 1 d 2 H{9) 1 d 2 Z (z) 



ri? (r) dr \ dr J r 2 H(6) d6 2 Z (z) dz 2 

where A is constant since the l.h.s. depends only on r, 9 while the middle depends 
only on z. 

The function Z (z) satisfies 

d 2 Z 



The BCs at z = 0, L imply 



+ XZ = 



= u(r,6,0) = R(r)H(0)Z(0) 
= u(r,9,L)=R(r)H(9)Z(L) 

To obtain non-trivial solutions, we must have 

Z (0) = = Z (L) . 

As we've shown many times before, the solution for Z (z) is, up to a multiplicative 
constant, 

/nnz\ /nn\ 2 
Z n (z) = sin J , A n =(^— J , ra = 1,2,3, ... 

Multiplying Eq. (92) by r 2 gives 

" (93) 



R (r) rfr V J H{9) d6 2 

where \i is constant since the l.h.s. depends only on r and the middle only on 9. 
We assume g (9, z) is smooth and 27r-periodic in 9. Hence 

H'(-tt) = H'(tt), H'(-tt) = H'(tt) 

We solved this problem above and found that /i = m 2 for m = 0, 1, 2, H (9) = const 
and 

H m {9) = A rn cosm9 + B m sin m6>, m — 1,2,3, ... 
Multiplying (93) by R (r) gives 
d ( dR(r 



r ^( r ^)"^ r2 + m2 ) i?(r) = 
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This is the Modified Bessel Equation with linearly independent solutions I m (ry/X^j 
and K m (r y/X^j called the modified Bessel functions of order m of the first and second 
kinds, respectively. Thus 



Since the J m 's are regular (bounded) at r = 0, while the K^s are singular (blow up), 
and since R mn (r) must be bounded, we must have c 2m = 0, or 

R (r)-r T (r,f\~\ - r T (—\ 

±L mn y ) °lm J m ^' V y n J °lm-'m y ^ y 

Thus the general solution is, by combining constants, 

oo oo 

u E (r, 9, z) = ^2 Im {~J~) sin ^ Amn C0S ( m6 ^ + Bm sin ( m9 ^ 



m=0 n=l 



In theory, we can now impose the condition u (a, 0,z) = g {9, z) and find A mn , B n 
using orthogonality of sin, cos. 
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